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Adrian Auer* and Guido Burkard 

Department of Physics, University of Konstanz, D-78457 Konstanz, Germany 

Defect centers in diamond are exceptional solid-state quantum systems that can have exceedingly 
long electron and nuclear spin coherence times. So far, single-qubit gates for the nitrogen nuclear 
spin, a two-qubit gate with a nitrogen-vacancy (NV) center electron spin, and entanglement between 
nearby nitrogen nuclear spins have been demonstrated. Here, we develop a scheme to implement a 
universal two-qubit gate between two distant nitrogen nuclear spins. Virtual excitation of an NV 
center that is embedded in an optical cavity can scatter a laser photon into the cavity mode; we 
show that this process depends on the nuclear spin state of the nitrogen atom. If two NV centers are 
simultaneously coupled to a common cavity mode and individually excited, virtual cavity photon 
exchange can mediate an effective interaction between the nuclear spin qubits, conditioned on the 
spin state of both nuclei, which implements a universal controlled-Z gate. We predict operation 
times below 100 nanoseconds, which is several orders of magnitude faster than the decoherence time 
of nuclear spin qubits in diamond. 


Introduction. Substantial experimental progress has 
been made in demonstrating the viability of nuclear spins 
coupled to nitrogen-vacancy (NV) centers in diamond as 
qubits. Compared with the NV electron spin, the nuclear 
spin offers for superior coherence properties, but so far, 
a scheme for the necessary two-qubit gates is lacking. 
Candidate nuclear spins are the intrinsic nitrogen nu¬ 
clear spin ( 14 N or 15 N) [1] or incidental proximal nuclear 
spins (e.g. 13 C) [2]. Decoherence times of TJ « 5 ms at 
room temperature have been measured [3], and elemen¬ 
tary single-qubit operations were implemented, includ¬ 
ing manipulation [1, 4-6], initialization [1, 4, 7-9] and 
high-fidelity single-shot readout [1, 8-11]. It was further 
demonstrated that the nitrogen nuclear spin can be a 
functioning part of a small quantum register [3, 9, 11 
14], or can act as a quantum memory to store and later 
retrieve the NV electron spin state [15]. Nuclear spin en¬ 
tanglement has been studied both experimentally [14,16- 
19] and theoretically [20-23]. However, a deterministic 
long-distance coupling scheme that does not utilize prior 
electron entanglement has not yet been demonstrated. 
The coupling of nuclear spins is fundamentally required 
in the context of quantum information processing, e.g. to 
perform universal quantum computation [24]. 

In this article, we develop and analyze a mechanism 
to optically generate a controlled quantum gate between 
two distant nitrogen nuclear spins (Fig. 1). The coupling 
between the nuclear spins is achieved by exchanging vir¬ 
tual cavity photons among two NV centers. External 
laser photons incident on each NV center can be scattered 
into the cavity mode, or vice versa, by exciting electronic 
Raman-type transitions between the NV ground and ex¬ 
cited state. We find that in the appropriate parameter 
regime, the scattering process depends on the nitrogen 
nuclear-spin state and can be completely suppressed for 
a specific nuclear spin configuration by properly tuning 
the laser frequency. This nuclear-spin dependent scat- 



FIG. 1. (a) NV center formed by a substitutional nitrogen 
(N) atom next to a vacancy (V) in the diamond lattice, (b) 
Magnetic field (B) direction along the defect axis (z) and 
N nuclear spin I. (c) Interaction between two NV centers 
(NV 1 and NV 2) that are coupled to the same mode io c of 
an optical cavity. The NV centers are excited by two lasers 
of frequency io L . Scattering of a laser photon mediates an 
effective interaction between NV 1 and NV 2. (d) Hyperfine 
levels mi = 0,±1 of the ms = —1 subspace for 14 N in the 
ground (gs) and excited (es) states, qubit states are indicated. 
S L and 5 C are the detunings of the laser frequency u> L from the 
mi = 0 orbital transition energy E s — A and from the cavity 
frequency io c . (e) Hyperfine levels mi = ±1/2 for 1S N. 
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tering mediates an effective interaction between two ni¬ 
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universal controlled-2T (cz) gate is implemented, which 
is equivalent to CNOT up to single qubit operations. A 
quantitative analysis of the proposed mechanism yields 
gate operation times below 100 nanoseconds, which is 
more than four orders of magnitude shorter than the de¬ 
coherence time of several milliseconds for the nitrogen 
nuclear spin. While cavity-mediated coupling between 
NV electron spins relies on the zero field splitting [25], 
the coupling of NV nuclear spins has its physical origin 
in the hyperfine interaction. 

Model. We start our analysis by describing a single NV 
center coupled to a single cavity mode and to an external 
laser field. The extension to two NV centers interacting 
with the same cavity mode, as required for the two-qubit 
gate, is straightforward and will be given later. To model 
the combined system of a single NV center, an optical 
cavity and the external laser, we use the time-dependent 
Hamiltonian 

H(t)=H NV + H c + H L (t ), (1) 

where H NV = H e +H n +Hh{ describes the electron (e) and 
nuclear (n) spin systems coupled through hyperfine (hf) 
interactions, H c the coupling to the cavity, and H L (t) 
the interaction with the laser field. In the presence of 
an external magnetic field B = Be z along the defect 
symmetry axis (z axis), the electron spin (S) and nuclear 
spin (I) Hamiltonians are given by (H = 1) [25-27] 

H e = le BS z + DSl - ^A S 2 z t z + 1 E g r z , (2) 

H n = - 7 n BI Z + QI 2 z . (3) 

Here, 7 e /27r = 2.803 MHz/G is the electron gyromag- 
netic ratio, E g = 1.945 eV is the energy gap between 
ground and excited state, and D — ( D gs + D es )/ 2 and 
A = D gs — D es with the zero-field spin splittings of 
the ground (D gs /27r = 2.88 GHz) and excited state 
(D es /27T = 1.42 GHz). The nuclear gyromagnetic ratio 
is denoted 7 „ and Q is the nuclear electric quadrupole 
coupling (see Tab. I). To describe the orbital degree 
of freedom, we use Pauli matrices r, ; (i = x,y,z), and 
choose the ground and excited states as t z eigenvectors 
with eigenvalues —1 and +1, respectively. We only con¬ 
sider the lower orbital branch of the excited state doublet 
( E y ). This is justified by naturally occurring strain fields 
of 10 GHz and more [28], which split the excited state 
into two well-separated orbital branches. 

Hyperfine interaction in the excited state is modeled by 
a diagonal hyperfine tensor, which has the same form as 
in the ground state [29, 30]. However, since the electron 
density at the nitrogen site is larger in the excited state 
[31], the hyperfine interaction is about 20 times stronger 
compared to the ground state according to measurements 
under ambient conditions [4, 32]. The difference 6A be¬ 
tween the hyperfine coupling in the ground and the ex¬ 
cited state forms the basis of the nuclear-spin dependent 
light scattering effect which we predict. Working at mag¬ 
netic field strengths away from the ground and excited 


TABLE I. Relevant nuclear-spin parameters for the NV cen¬ 
ter. 


Parameter 

14 n 

1S N 

Nuclear spin I 

1 

1/2 

7n/27T 

0.308 kHz/G [33] 

-0.432 kHz/G [33] 

Q/2tt 

-5 MHz [1, 4, 34] 

0 

A ss /2n 

-2.2 MHz [1, 4, 5, 34] 

3.0 MHz [32, 34] 

A es /27T 

40 MHz [4] 

61 MHz [32] 


state level anticrossings, electron-nuclear spin flip-flop 
processes are energetically suppressed. Therefore, we ne¬ 
glect the transverse part of the hyperfine tensor and only 
include the longitudinal coupling. Denoting the hyper¬ 
fine coupling strengths by A gs and A es for the ground 
and excited state (Tab. I), we arrive at 

H h( = AS Z I Z + ^ 6A t z S z I z , (4) 

where A = (A es + A gs )/2 and 5 A = A es — A gs . 

We consider the NV electronic orbital transition be¬ 
tween the ground and excited state to be coupled to a sin¬ 
gle mode of the optical cavity, which, in the rotating-wave 
approximation, is described by H c = co c a'a + g(r+a + 
r_a^), where co c is the cavity frequency, is cavity- 
photon annihilation (creation) operator, g the coupling 
strength (which be assumed real), and t± = (T x ±ir y )/ 2. 
The external laser is described by a classical field of fre¬ 
quency w L that excites electronic orbital transitions be¬ 
tween states having the same spin projections mg and 
to/, H h (t) = fle~ lu>Lt T- (_ + H*e lWLt r_. Here, Q is the 
complex Rabi frequency that depends on the phase of 
the laser field. The Hamiltonian H(t) can be made 
time-independent by transforming into a rotating frame, 
H' = e&Htye-** - £ with £ = uj h (a^a + r z /2), and we 
obtain H' = H' e + H n + + H' C + H[. The transformed 

part H' e of the electronic Hamiltonian is obtained by re¬ 
placing E g with the detuning 6 h = E g — w L in H e . In the 
Hamiltonian H c , the transformation causes a shift of the 
cavity frequency to S c = w c — w L , which is the detuning of 
the laser from the cavity mode. The laser Hamiltonian 
H L (t) becomes time-independent, H[ = Qt + + H*r_. 

Nuclear-spin dependent scattering. Virtually exciting 
the NV center by the external laser field can finally lead 
to an excitation of the cavity mode through the coupling 
g. We describe this process by using quasi-degenerate 
perturbation theory in terms of a Sclirieffer-Wolff (SW) 
transformation [35, 36] to eliminate the intermediate vir¬ 
tual transition to the excited state, and obtain a model 
that effectively describes the scattering of a laser photon 
into the cavity mode, and vice versa, that particularly 
depends on the nitrogen nuclear spin projection ?n/. It 
is exactly this spin-dependent scattering that eventually 
enables a conditional two-qubit quantum gate. 

To implement the SW transformation, we construct 
an anti-Hermitian operator S such that [5, H 0 ] = V (see 
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Appendix), where the part Hq = H' e + H n + H^i + <5 c aJa 
only acts on the ground and excited state manifold, re¬ 
spectively, and V describes transitions between these 
two Hilbert subspaces. In the transformed Hamiltonian 
H = e~ s H'e s , we keep the lowest order in the inter¬ 
action V and continue with the effective Hamiltonian 
H « H 0 + [V,S\/2. The effective ground-state Hamil¬ 
tonian becomes (see Appendix) 


(gB) = - (^gs + 7n B) I z + Qll + <S c ot 

1 


2 L' 


( 5 ) 

g n ({5 A I z - (5 l ) -1 + (6A I z + S c - 6J- 1 ) a* + h.c. 


Here, we restrict our consideration to the ms = — 1 sub¬ 
space, and define the detuning d L = d L — A of the laser 
frequency from the mj = 0 orbital transition (Fig. 1). We 
omit all constant terms and neglect small energy shifts 
proportional to g 2 and |S2| 2 . 

On the basis of previous experimental work [5, 14, 15] 
using the 14 N nuclear spin as a qubit, we choose the nu¬ 
clear spin sublevels |m/ = 1) = |1) and |mj = 0) = |0) as 
the computational basis. We can neglect the mj = — 1 
state because the transition frequency between these two 
levels is well separated from other transitions [5]. From 
Eq. (A8), one can see that the effective coupling of the 
NV center to the cavity via the virtual laser excitation 
depends on the spin state of nitrogen nucleus and can, 
e.g., be completely suppressed for one of the two spin 
states. This is the case if the laser frequency is cho¬ 
sen such that e.g. S L = S c /2 , where only scattering from 
the mi = 1 state is possible. By using I z = 11) (11 and 
1 = 11)(11 + |0)(0|, we find the qubit Hamiltonian 

ff (qubi t) = (Q- A gs - ln B) |1)<1| + 6 c a*a 

+ <7 , |l)(l|a t + (<7 , )*|l)(l|a, (6) 


with an effective coupling strength 


g' = gil 


5A 

^ 2 -d) 2 ' 


( 7 ) 


Scattering only from the mi = 0 state is possible for 
<5 l = 5A + 5 C / 2 occurring with the same effective coupling 
strength g' [Eq. (7)]; however, we concentrate on mi = 1 
scattering in the following. 

Spin-spin interactions. To understand the scattering 
mechanism of a laser photon into the cavity mode qual¬ 
itatively, we so far neglected spin-mixing terms in the 
lower branch of the excited state doublet [26, 32, 37 
39]. However, to make quantitative predictions of the 
effective scattering process, we take into account the fine 
structure of the excited state manifold. So far, electronic 
spin-spin interactions were only incorporated by the zero- 
field splittings D gs . and D es . In the limit of high strain 
considered here, the two branches of the excited-state or¬ 
bital doublet split and anticrossings in the lower branch 
mix spin states with different quantum numbers ms- The 



FIG. 2. Ratio /(<5 C ) = gig' of coupling strengths with (g) 
and without (</) spin-spin interaction in the excited state. 
Magnetic field strengths B are chosen such that mi is a good 
quantum number in the ground and excited state. 


Hamiltonian describing the spin mixing is [27, 39] 

H s = \(1 + r z ) {Si - S 2 y ) - ^ {S X S Z + S z S x )j , 

( 8 ) 

where transitions between the excited state orbitals have 
been neglected due to the high strain, and the fine struc¬ 
ture parameters are given by Ai/27t = 1.54 GHz and 
A 2 / 27 r = 0.154 GHz [39]. The effective coupling strength 
g analogous to Eq. (7) can be obtained by adding H s to 
the bare NV Hamiltonian H NV , and then performing the 
SW transformation. In doing so, we assume the cavity 
to be populated by at most one photon, and only if the 
NV center is in the ground state. In the excited state, we 
need to include all spin states ms = 0, ±1. The effective 
ground state Hamiltonian in the case of to/ = 1 scatter¬ 
ing has the same form as given in Eq. (6) with a different 
coupling strength g = g'f(S c ). The detuning-dependent 
part f{S c ) is plotted in Fig. 2. 

Controlled quantum gate. For the two-qubit gate, we 
consider two NV centers (z = 1,2) coupled to the same 
cavity mode and each individually driven by a laser of 
frequency w L [Fig. 1 (c)]. In the following, we keep only 
the lowest order of the interaction parts, and consider 
m/ = 1 scattering on both NV centers. Furthermore, 
we assume detunings d L and S c such that the cavity is 
excited only virtually, which, in turn, leads to an ef¬ 
fective interaction between the two NV centers. To de¬ 
scribe this interaction, we apply a second SW transforma¬ 
tion to H ( 2 gs) = 5 c o t a + X) 2 =1 (Q - Ag S - j n B) |l)j(l| + 
(gi|l)i(l|a 1 '+h.c.) to eliminate the cavity mode by choos¬ 
ing S = — X^=i(ff*/ <! k|l)i(l|a' i ’ — h.c.), which leads to an 

effective Hamiltonian H c s = e~ s H 2 S ^ e s , where again 
only the lowest order contribution of the off-diagonal el¬ 
ements is kept (see Appendix). H e s comprises single¬ 
qubit terms = {Q - A gs -"/ n B - |g i | 2 /<5 c )|l) i (l|, and 
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FIG. 3. (a) Effective two-qubit coupling strength g\ 2 as 

a function of 5 C for different values of the laser Rabi fre¬ 
quency fii = = D (see legend, valid for all figures) and 

gi/2ir = gi/2-K = 50 MHz for 14 N nuclear spins, (b) Time 
Tcz to generate a CZ gate between the two nuclear spins as 
a function of 5 C for gi/2n = g 2 /2n = 50 MHz. (c) and (d) 
equivalent for 15 N nuclear spins. All calculations performed 
at B = 120 G. 


a two-qubit interaction term, 

#i„t = -<?i2|ll>(ll|. (9) 

Here, 111) = 11) 1 1 1)2 is the nuclear spin state of both 
NV centers 1 and 2, and the effective two-qubit coupling 
strength gi2 is found to be 

ff!2 = 2 ^ I C0S(^1 — (f> 2 ), (10) 

where </,■ denotes the phase of the ith laser field, f1, = 
. Quantitative predictions of g 12 are plotted in 
Fig. 3 (a). 

Since [H^g,H- lnt ] = 0, the time evolution U generated 
by the Hamiltonian H e s can be written as 

U(t) = e~ iHeCit = (U^t) ® U 2 (t)) U 12 (t), (11) 

where 17* (i) is a single-qubit rotation of nuclear spin i and 
U\ 2 (t) describes a two-qubit operation generated by the 
interaction part H lnt . In Eq. (11), the time evolution of 
the cavity field has been omitted, since the nuclear spin 
degree of freedom has been decoupled from the cavity 
field by the above transformation. In the following, we 
only concentrate on the two-qubit interaction part, and 
disregard single-qubit rotations since they can be undone 
afterwards, e.g., by off-resonant excitation of the ground- 
state electronic spin transition, thereby implementing a 
phase gate on the N nuclear spin [5], or direct driving of 
the nuclear spin transitions [40]. 

For an operation time of r cz = Tf/gi 2 , a CZ gate is 
implemented on the two nuclear spin qubits, 

U 12 (t cz ) = |00)(00| + |01)(01| + |10)(10| - |11)(11|, (12) 


from which CNOT can be created using additional 
Hadamard gates [24], In Fig. 3 (b), values of r oz are 
shown for different Rabi frequencies Q. As the main re¬ 
sult of our paper, we find fast operation times clearly 
below 100 ns. In our calculations, we assumed large de¬ 
tunings |<5 l — <5 C | > g and |<5 L | > |H| to justify the effective 
model used. 

Conclusions. Nitrogen nuclear spins in diamond have 
proved to be highly promising candidates to physically 
realize qubits. We have presented a theoretical proposal 
for the implementation of a controlled optical cavity- 
mediated quantum gate between two nitrogen nuclear 
spin qubits intrinsic to NV centers in diamond. Gate op¬ 
eration can be achieved within 100 nanoseconds or less, 
which is more than four orders of magnitude below the 
nuclear-spin decoherence time. Assuming r cz ss 20 ns [see 
Fig. 3 (b)], the cavity loss rate k must not exceed values 
of 1/t cz ~ 50 MHz, which requires Q factors of 10 6 10' 
[25]. High-Q silica microsphere cavities can reach such 
values [41] and progress towards photonic crystal cavi¬ 
ties in bulk diamond exceeding Q-factors of 10 5 has been 
recently achieved [42]. 

In addition to the presented findings, an equivalent 
analysis for the 15 N nuclear spin with I = 1/2 show that 
the proposed scheme also works for this isotope if the 
computational basis is chosen as |1) = | mj = +1/2) and 
|0) = |to/ = —1/2) [Fig. 1 (e)]. We find the same effec¬ 
tive scattering rate g' [Eq. (7)] for mi = ±1/2 scattering 
for laser detunings <5 L = (<5 C ± 5A)/2. Including spin-spin 
interactions, the effective two-qubit coupling strength g V2 
and the gate time r cz show qualitatively the same behav¬ 
ior as for the 14 N nuclear spin, and are depicted in Figs. 3 
(c) and (d). During the fast electronic excitation cycles, 
the nuclear spins are subject to a time-varying hyperfine 
interaction. However, using a spin-fluctuator model, it 
has been shown that nuclear spin state will be unaffected 
and coherence can be preserved [43]. Together with el¬ 
ementary and experimentally demonstrated single-qubit 
operations, the realization of a universal CZ gate makes 
the nitrogen nuclear spin valuable for quantum compu¬ 
tation in addition to its remarkable quality as a quantum 
memory [15]. 
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Appendix A: Schrieffer-Wolff transformation to 
eliminate excited state 

We separate the Hamiltonian H' into a block-diagonal 
part Hq that only acts within the ground and the excited 
state manifold, respectively, and an off-diagonal part V 
that connects these two manifolds, 


H'= H 0 + V. (Al) 

To implement the Schrieffer-Wolff (SW) transformation 
[35, 36], we construct a unitary transformation exp(— S) 
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with some anti-Hermitian S to obtain a new Hamiltonian 

H, 

H = e~ s H'e s , (A2) 

that contains no matrix elements that connect the ground 
and the excited states up to a desired order in V. If we 
choose the anti-Hermitian operator S in such a way that 

[S,H 0 ] = V (A3) 

holds, the leading order in V cancels. If we keep the 

lowest order in V, the Hamiltonian H is approximately 
given by 

H**H Q + hy,S\. (A4) 

The block-diagonal part Hq of H' is given by 

H 0 = H' e + H n + H h( + <5 c a t a, (A5) 

and the interaction terms are 

V = g{r+a + r_a^) + f2r+ + fl*r_. (A6) 

From the condition in Eq. (A3), we find 

s = n (a S 2 - A h{ S z I z - ~s L ) _1 T+ - h.c. 

+ g [asI - A.mS z I z + 5 C - T+a- h.c., (A7) 

and the effective Hamiltonian for the decoupled ground 
state manifold becomes 

H (gs) = - (A gs + ln B) I z + QI 2 + S c a^a 

+ i [ g n ((A hf I z - <> L ) -1 + (A hf I z + S C - 50” 1 ) ^ 

+ h.c. . (A 8 ) 

Here, we restrict our consideration to the ms = — 1 sub¬ 
space, and define the detuning S L = S L — A of the laser 
frequency from the mi = 0 orbital transition. We omit 
all constant terms and neglect small energy shifts pro¬ 
portional to g 2 and |fl| 2 . 

Appendix B: SW transformation to eliminate virtual 
photon 

We start from a Hamiltonian H^ that describes two 
NY centers (* = 1,2) coupled to a common cavity mode 
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